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Abstract 

An analysis of pseudorapidity correlations in 10.7 A Gev 197 Au and 200 A 
GeV 32 S interactions with emulsion nuclei by the normalized range method has 
been carried out. The evidence for events with large multiparticle correlations 
is presented. The most significant "correlation force" effect is apparent in the 
interactions of light emulsion nuclei (CNO group) and gold nuclei at an energy 
of 10.7 A GeV, corresponding to an absolute disintegration of the target nucleus 
n b + n 3 ~ 0. 

1 Introduction 

A study of interactions of relativistic nuclei is related to a search for unusual nuclear 
matter conditions (as quark-gluon plasma (QGP), for example), which can be observed 
at high energies [ 1 1. The behaviour of simple features of these processes, such as to- 
tal cross sections, multiplicity and pseudorapidity distributions, etc. is rather precisely 
described by various models, even if they are based on the various assumptions about 
nuclear structure and dynamics of nuclear interactions. More detailed information, crit- 
ical to a choice between existing models, can be deduced by investigating a fluctuation 
structure of secondary particle distributions. 

Recent studies of pseudorapidity distributions have shown that some considerable 
deviations from average distribution can be observed in individual events. Certainly, 
these deviations can be initiated by the statistical factors, connected with finite number 
of particles in the given event. If the statistical explanation cannot describe these fluc- 
tuations, there are the some dynamic mechanisms, which could produce peculiarities 
of pseudorapidity distributions. 

In our previous work [2 1 we have proposed a normalized range method for analysing 
pseudorapidity correlations in multiparticle production processes. This method allows 
not only to discriminate dynamic correlations from statistical ones, but also to deter- 
mine a "force" and "length" of these correlations. 

In the present paper we apply this method to the analysis of experimental pseudo- 
rapidity distributions of secondary particles obtained in 10.7 A GeV 197 Au and 200 A 
GeV 32 S interactions with emulsion nuclei. 



2 Analysis procedure 

As result of high-energy interaction of two nuclei plenty of secondary particles is pro- 
duced. According to the existing notions, secondary particles, which are "emited" 
from "interaction volume", have pseudorapidities, corresponding to a central region 
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of pseudorapidity distribution. At borders of the distribution fragments (of the target- 
nucleus and projectile-nucleus) bring in considerable contribution. And so, in order to 
investigate of pseudorapidity correlations in distribution of particles from "interaction 
volume" we have chosen pseudorapidity interval A77 = 4 (for the 32 S, 200 A GeV, 
Vmin = 1-0, Vmax = 5.0 and for the 197 Au, 10.7 A GeV, r\ min = 0.3, rj max = 4.3). 
This interval has been subdivided into k parts with Si] = Arj/k, By counting the 
number of particles in each subinterval we arrive at a sequence n\. A pseudorapidity 
fluctuation, or the normalized relative deviation of an individual event from average 
pseudorapidity distribution Q is given by 

= n e Jn e -nj/n = rif n_ _ 1 ^ 
Ui/n n e rii 

where n\ is the number of particles in the i-th bin of an event with particles number 
n e , and n< = ^ e nf is the total number of particles for all events in the i-th bin, and 
n = J2 e nC i s tne tota l num ber of particles for all events. 

For analysis we have selected events with the number of secondary particles (ex- 
cept fragments) greater than 90, because high statistics are necessary for a study of 
correlations by our method. 

As described in our previous work |]2] for an investigation of pseudorapidity cor- 
relations we analysed the normalized range H(k') — R(k')/ S(k') (where R(k') and 
S(k') are a "range" and a standard deviation, which is calculated by Eqs.©-©, see 
below) versus the size of the pseudorapidity interval dr/ = k'Srj, (1 < k' < k) using a 
function 

H(k') = (ak'f (2) 

where a and h are two free parameters and h is the correlation index (or Hurst index). 
If the signal & represents white noise (a completely uncorrelated signal), then h = 0.5. 
If h > 0.5, long-range correlations are in a system 

In our calculations we have used k = 8192. The choice of such large value of k is 
not necessary (it is possible using lesser k also). But the more k the more accuracy of 
method is approachable. 

So, for the sequence £j, 1 < i < k, quantities of R(k) and S(k) were calculated by 
following formulas: 

1/2 
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R(k) — maxX(m, k) — minX(m, k) (4) 

l<m<k l<m<k 

where the quantity X(m, k) characterizes the accumulated deviation from the average 

k 



^ k 



'it is possible to use also absolute deviation f j = n^/n e — rii /n. 
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for a certain interval mSr/, 

m 

X(m,k) = ^[&- < f >], l<i<m<k (6) 

z=l 

Then the sequence £j has been subdivided into two parts: ^, 1 < i < k' = k/2 
and k' + 1 < i < k, and the value of #(fc/2) = R{k/2)/ S{k/2) was found 
for each of the two independent series. Similarly ^ and have been subdivided 
further, followed by the calculation of H(k/A) = R(k/4)/ S(k/4). This subdivision 
and analysis procedure for newly obtained series of ^-values is continued until dr\ > 
0.1. H corresponding to the same value of k 1 have been averaged and drawn on a 
log-log scale as a function of k' (Fig.l). 




ln(k') 



Figure 1: The correlation curve for (Au+Em)-interactions at 10.7 A GeV. Points are 
calculated values, lines are fits by © 

For our correlation analysis it is possible to use this dependence visually and de- 
termine the points belonging to the same line and their fit by (O. But more precisely 
and more clearly one can investigate the behaviour of the correlation index h versus 
the width of pseudorapidity bin drj with the help of a 3-point fitting procedure. Those 
three points were: the first has been subsequently conferred by obtained Hurst index 
and other two ones are neighboring points. For example, h corresponding to k' = 1024 
(In(k') = 6.931, di] = 0.5) in the Fig.l has been found, using and with following 
a fitting procedure for points k' = 512 (In(k') = 6.238), k' = 1024 (In(k') = 6.931), 
k = 2048 (ln(k') = 7.625) from the Fig.l. 

3 Data 

The above analysis procedure has been applied to three samples of events 
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1. Experimental data of the 32 S 200 A GeV (264 events) and the 197 Au 10.7 A 
GeV (315 events) interactions with emulsion nuclei. NIKFI BR-2 stacks of nuclear 
emulsions, sizes of (10 x 20 x 0.06) cm 3 have been exposed horizontally to the 10.7 A 
GeV 197 Au beam at the BNL/AGS and the 200 A GeV 32 S at the CERN/SPS. Details 
about scanning and measurements are given elsewhere ||5l,||6l. 

2. "Data" obtained by FRITIOF 7.02 model Q, of the 197 Au (10.7 A GeV, 315 
"events") and 32 S (200 A GeV, 264 "events") interactions with emulsion nuclei. 

3. A random process simulated by the HRNDM program, included in the program- 
ming package HBOOK |8]. Multiplicities and pseudorapidity distributions for this 
process are similar to the experimental data (Au + Em at 10.7 A GeV, 315 "events" 
and S + Em at 200 A GeV, 264 "events"). 



4 Results and discussion 

It should be noted that the £j distributions for all data samples are similar (at least 
formally) to themselves. At the same time, the experimentally observed dependence of 
the correlation index on the size of the interval dr] is different from both HRNDM and 
FRITIOF events (see Fig. 2). 



■ Experimental data 
A FRITIOF data 
* Stochastic data 
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Figure 2: The value of Hurst index h versus a width of pseudorapidity bin dr) for data 
of Au+Em 1 0.7 A GeV (left) and S+Em 200 A GeV (right) 



As it is seems from Fig. 2 the most obvious difference is evident at large dr) for 
(Au+Em) interactions at 10.7 A GeV, where a significant increase of the correlation 
index is observed, thus indicating stronger correlation "force" Q. 

A detailed analysis of individual events has shown that all experimental events can 
be conditionally divided into two classes. In Fig. 3 correlation curves for two events 
from two different classes, is presented. 

Firstly, the events without a "kink" of the correlation curve correspond to random 
pseudorapidity distribution of secondary particles. In these events fluctuations are es- 
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Figure 3: The correlation curve for two of events of (Au+Em)-interacions at 10.7 A 
GeV 



sentially statistical, connected with the finite number of particles in the event. Secondly, 
the events "with kink" of the InH(k') versus Ink' correlation curve, in which we do 
observe significant multiparticle correlations at Srj > 1. 

The event selection was made by the average correlation index 

„„ = » '£ WWl-W*-,)) , lki _ km (7) 

imax - 1 in(ki) - ln{ki-i) 

which was determined by a power-law fit (fJJ to the correlation curve curve InH(k') 
versus Ink'. If h ev was greater than 0.62, then the event was related to the group of 
events with kink. We can analyze the more correlated events by increasing the criterion 
of selection up to h ev > 0.7 or h ev > 0.8. In Fig.4 the curves of InH(k') versus Ink' 
for 0.62 criterion are shown. 

The behaviour of these without kink curves are similar to the FRITIOF events. 

It is interesting to note that the relative number of correlated (S+Em)-events was 
considerably lower than for (Au+Em) interactions. Namely, we have 187 events with 
h ev > 0.62 for gold and only 97 ones for sulfur. The (S+Em)-interactions with very 
strong correlations (with h ev > 0.8) were absent at all. In the case of (Au+Em)- 
collisions the number of such events was 16. 

To study possible distinctions in the mechanism of formation of a final condition 
for the two types of interactions (with kink and without kink of the correlation curve), 
we have analysed the behaviour of fragments of the projectile and the target nuclei. 
Thus, most difference is revealed for the Nh -distributions, /i-particles are called a 
sum of b- and g-particles (Nh = n& + n g ). According to the existing photoemulsion 
experiments criteria the fe-particles are fragments of target nucleus with kinetic energy 
E kin < 26AM eV, g-particles are recoil protons with 26 < E kin < AOOMeV. @ As 



2 It should be noted that h-particles and Hurst index (h) have not any connection between themselves. The 
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Figure 4: The correlation curve for two of events of (Au+Em)-interacions at 10.7 A 
GeV The value of Hurst index h versus a width of pseudorapidity bin drj for the events 
with h e v > 0.62 and the events with h ev < 0.62 from experimental data of Au+Em 
10.7 A GeV (left) and S+Em 200 A GeV (right). 



it is seen from Fig. 5 a large part of events with kink, as opposed to the interactions 
without kink, proceed with a complete (or almost complete) disintegration of the target 
nucleus (peak in area Nh = 0). 

Most of these kink events (events with very strong correlations) corresponds to such 
events (Fig. 61eft). 

The photoemulsion consist of Hydrogen(39.2%), nuclei of CiVO-group(35.3%) 
and AgBr nuclei (25.5%). In order to know what interactions give main contribution 
in peak of Nh = we have analysed n s distribution for events with rib = (Fig. 
6right). 

The first peak in this distribution at n s < 30 corresponds to peripheral interactions 
of gold with photoemulsion nuclei and central interactions, (p + Au) . The second peak 
at 30 < n s < 200, is a result of central (since n s are large) interactions of light nuclei 
(CNO group) and nuclei of gold. The peak in the area of n s > 200 corresponds to 
central collisions of (Ag + Au) and (Br + Au). 

Hence, by analysing events with n s > 200, we study interactions of heavy nuclei 
of Ag and Br with nuclei of gold. We expect that greatest correlations in such in- 
teractions. However, as it can be seen from Fig. 71eft, the peak in area of Nh = 0, 
corresponding to the greatest values of a correlation index, has disappeared! 

Therefore, the events which give the main contribution to the peak in the area of 
Nh ~ and which have great kink of correlation curve, correspond to central interac- 
tions of Au and light nuclei of CNO group. 

It is possible that there are several probable reasons for the existence of such large 

signs have coincided by accident. 

3 s-particles are secondary ones with E kin > 400AM eV . 
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Figure 5: The iV^-distribution for data of Au+Em (10.7 A GeV) for the events with 
h ev > 0.62 (left) and for the events with h ev < 0.62 (right) 



correlations in multiparticle production processes. But, nevertheless, at present we 
know only one way when similar behaviour (to experimental one) of the correlation 
curve can be obtained. As result of interactions of two particles at high energy there is 
formed some large claster (or several clasters), which then disintegrates into a plenty 
of secondary particles in the interval of Sr) < 1. 

This process was simulated with the help of the programme "HRNDM", included 
into the software package of "HBOOK" [8]. We simulated stochastic events with pseu- 
dorapidity distributions similar to experimental data. We have designate: a total num- 
ber of generated "particles" in an event as n e , probability of production of claster as p 
and number of the clasters as j (0 < pj < 1). Each simulated "particle" can be dis- 
integrated (with probability pj) into pn e -particles with pseudorapidity rji = rj m + fi, 
1 < I < pn e , where fi is random function, which changes with uniform probability in 
region [-RN,+RN], r\ m is pseudorapidity of the "parent particle". 

By varying p, j and RN, we can define curves most precisely reproducing the 
behaviour of our experimental data. 

From Fig.7right it is seen that the more dimension of the compound system, the 
more "kink" of the correlation curve will be displayed. 

Moreover, such the large kink of correlation curve, which observed for experimen- 
tal events, cannot be explaned by two-particle disintegrations. So, the curve with j=50 
and p=0.02 in Fig.7right corresponds to process with maximal two-particle correlations 
(all secondary particles are produced as result of two-particles disintegrations). In this 
case the correlation index reached only 0.62. In order to obtain the values h which 
observed for experimental data, it is necessary to assume that the compaund system, 
which produces ~ 40 % (for Au+Em) and ~ 24% (for S+Em) of secondary particles, 
is formed. 

It should be noted that for search of the most interesting effects as quark-gluon 
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Figure 6: The iV^-distribution for data of Au+Em (10.7 A GeV) for the events with 
h ev > 0.8 (left) and the iV s -distribution for data of Au+Em (10.7 A GeV) for the 
events with rib — (right) 



plasma, for example, in CERN experiments as NA44, NA49, WA98 etc. symmetric 
collisions of heavy ion as Pb+Pb (or Au+Au) are used. If our results will be confirmed 
by other experimental data (for instance, in interactions of Au (or Pb) and nuclei of 
CNO group at other (more or less than 10.7 a GeV) energies) then it will be possible 
declare that probability of exhibition of QGP in asymmetric interaction is more than 
that in symmetric ones. 

To clear this point it is possible by primitive modelQ. Supposing that a nucleon is an 
elastic ball, inside which quarks are disposed, and nucleus is some system consisting 
the balls, which are glued together. At the moment of collision first balls, coming 
into contact, are squeezed by other balls, which push them from behind and in front. 
And so, a real energy of interaction (at this "region") is significantly more than an 
energy per nucleon (some cumulative effect). A pressure (at this "region") is enormous. 
Volumes of the "interacting balls" is getting significantly diminished (is squeezed). 
As result a nucleon "envelope" can burst and quarks (from the nucleon) will be in 
quasifree state inside dense "encirclement" of squeezed (but not bursting) neighbouring 
nucleons. This region can exist and be expanded till large nucleus are not split (when 
"force of glue" will be less than "breaking force of the wedged (smaller) nucleus"). If 
an interaction is less central (almost peripheral) one, then quarks region and time of 
its existence will be less, because of that the probability to chop off smaller "scrap" is 
more than that to split nucleus to "fifty-fifty". In case of collision of symmetric nuclei 
the quarks region can appear also, but its dimensions will be still less, because dense 
"encirclement", which is capable to retain (inside itself) quarks, is absent (it is similar 

4 It should be noted that the process of nuclei interaction at high energies is more complicated one. Nev- 
ertheless, the model allow to realize why largest correlations should be exhibited in asymmetric collisions 
instead of symmetric ones". 
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Figure 7: The iV^-distribution for data of Au+Em (10.7 A GeV) for the events with 
n s > 200 (left) and the value of Hurst index h versus a width of pseudorapidity bin Ctq 
for modelled data with R N = 1 and different p and j . 



to peripheral interaction of asymmetric nuclei). 

5 Conclusion 

The analysis of pseudorapidity correlations in 10.7 A Gev 197 Au and 200 A GeV 32 S 
interactions with emulsion nuclei has been carried out by the normalized range method. 
This analysis has revealed events with large multiparticle correlations (with large kink 
of correlation curve). We assume that the main cause of this behaviour of the corre- 
lation curve is a formation of large compound system, which then disintegrates into a 
large number of secondary particles in the interval of 77 ± Srj (5r) < 1). The most sig- 
nificant pseudorapidity correlations are observed in the central interactions of nuclei, 
which have considerable difference in volume (in nuclear weight, charge, etc.) at rather 
"low" energies (the nuclei of Au and CNO group at 10.7 A GeV). 

It would be interesting to see whether a similar effect is observed in interactions of 
other nuclei (for example, lead and nuclei of CNO group) and whether it depends on 
the interaction energy (if the effect will be increase at energy of collision more than 
10.7 A GeV). 
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